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Abstract 

We study the quark mass dependence of various response functions, which 
contribute to chiral susceptibilities and the specific heat in the staggered fermion 
formulation of two-flavour QCD. This yields information about the critical expo- 
nents a, (3 and d. In the case of the chiral susceptibility, obtained as derivative 
of the chiral order parameter with respect to the quark mass, we calculate all 
contributions. This allows to construct a cumulant of the order parameter, which 
is a scaling function and yields a direct determination of the critical exponent 5. 
All our results are consistent with a second order phase transition. 



1) Introduction 



In contrast to QCD with three or more Ught flavour degrees of freedom, where nu- 
merical simulations seem to indicate that the finite temperature chiral symmetry restoring 
phase transition is first order, there is increasing evidence that the two-flavour theory has 
a second order phase transition in the limit of vanishing quark masses. At least in the 
parameter range studied so far in the staggered formulation of lattice QCD, no signs of a 
flrst order transition have been observed [1]. In fact, in a recent analysis [2] of the exist- 
ing Monte Carlo data for pseudo-critical couplings, gdrriq), at non-vanishing quark mass, 
one of us has shown that the functional dependence on niq is consistent with the scaling 
behaviour expected for a chiral phase transition with critical exponents in the universality 
class of a 3-d 0(4) symmetric cr-model [3,4]. Still, even in the class of 0(4) symmetric 
models the occurrence of first order phase transitions can not be ruled out easily [3,5]. A 
careful analysis of the nature of the chiral transition is thus asked for. 

In most of the existing studies of the chiral phase transition in the staggered fermion 
formulation of two-fiavour QCD [6,7] the critical parameters have been determined through 
a subjective inspection of the dependence of the chiral order parameter or the Polyakov loop 
expectation value on the gauge coupling. In Ref. [8] the volume and quark mass dependence 
of the critical couplings have been examined by locating the maximum of the variance of 
these quantities. They are, however, not directly related to bulk thermodynamic quantities 
and predictions for their scaling behaviour can, therefore, not directly be obtained from 
properties of the two-flavour partition function, ie. the singular part of the corresponding 
free energy. A priori it thus is not clear in how far these quantities can give information 
on the critical exponents of two-flavour QCD. We will discuss this question here in detail. 

In this paper we will study susceptibilities and the speciflc heat obtained from second 
derivatives of the partition function, 
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The determination of the maxima of these quantities will allow us to give a precise quanti- 
tative definition of the pseudo-critical couplings at non-vanishing values of the quark mass. 
We will discuss the scaling behaviour of these couplings as well as that of the peak heights 
in the response functions. 



From Xm cind the chiral order parameter 
we will construct the cumulant 

which allows a direct investigation of the zero-mass critical coupling as well as the critical 
exponent 6 [9]. In such an analysis it is particularly important to determine the complete 
chiral susceptibility, Xm, ie- taking into account also the contribution from the connected 
part in the product of the quark bilinears. 

This paper is organized as follows. In the next section we will summarize basic re- 
lations for the critical behaviour of thermodynamic quantities in the vicinity of a second 
order phase transition, which is controlled by the temperature as well as an external sym- 
metry breaking field, ie. the quark mass. In section 3 we discuss the observables from 
which critical exponents for QCD with dynamical quarks on the lattice are calculated. In 
section 4 our Monte Carlo data for various response functions is presented for three values 
of the quark mass on lattices of size 8^ x 4. The determination of the critical exponents 
a, (5 and 5 is discussed. Finally we give our Conclusions in section 5. 



2) Scaling Relations 



If the chiral phase transition in two- flavour QCD is second order, the critical behaviour 
of thermodynamic quantities will be controlled by two external parameters, the reduced 
temperature t = {T — Tc)/Tc and the symmetry breaking field h = niq/T. In the staggered 
formulation of lattice regularized QCD these two dimensionless couplings are given by 

^_ 6 6 



gm ' (4) 
h = rUgNj. . 

Here (/^(O) denotes the critical coupling on a lattice of temporal extent Nt- in the limit of 
vanishing bare quark mass. For non-vanishing values of the quark mass, a pseudo-critical 
coupling, g1{mq), can, for instance, be defined as the location of a peak in one of the 
susceptibilities or the specific heat defined in eq. (1). 

In the vicinity of the critical point the behaviour of bulk thermodynamic quantities 
is governed by thermal {yt) and magnetic [yh) critical exponents, which characterize the 
scaling behaviour of the singular part of the free energy density, 

f{t,h) = -^\nz = b-'f{by%by'^h) . (5) 



Here b is an arbitrary scale factor. It is expected that the chiral phase transition in QCD 
can be described by an effective, three dimensional theory for the chiral order parameter, 
which in the case of two-flavour QCD would amount to an 0(4) symmetric spin model 
[3,4]. Still the generic structure of this effective theory leaves open the possibility of a first 
order transition [3,5]. A quantitative analysis of the scaling behaviour of thermodynamic 
quantities is needed to further support the existence of a second order phase transition 
in the zero quark mass limit. In the staggered lattice discretization of QCD the situation 
gets even more involved due to the fact that flavour symmetry is partially broken by 0{a) 
terms, where a is the lattice spacing. Correspondingly, at least at large lattice spacings, 
the relevant symmetry group is 0(2) [9], which would lead to somewhat different critical 
exponents and even might influence the order of the phase transition [9] . It remains to be 
seen if and when the full flavour symmetry is effectively restored on lattices with small but 
finite spacings. 

The scaling behaviour of the specific heat and susceptibilities is controlled by the 
critical exponents a = {2yt — l)/yt, P = ~ yh)/yt and S = yh/{l — yh)- Their numerical 
values for 0(2) and 0(4) symmetric spin models in three dimensions are given in Table I. 
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0(2) 


-0.007(6) 


0.3455(20) 


4.808(7) 


-0.004(3) 


0.792(1) 


0.394(2) 


0(4) 


-0.19(6) 


0.38(1) 


4.82(5) 


-0.10(3) 


0.793(3) 


0.34(1) 



Table I: Critical exponents of 3-d 0{N) symmetric spin models [4,10,11]. 
The exponents Zm and Zt characterize the scaling behaviour of susceptibilities 
as defined in eq. (6). The exponent Za is defined in eq. (10) and describes the 
leading mass-dependence of the specific heat peak. 

Eq. (5) can be used to extract scaling laws for various quantities, valid in the vicinity 
of the critical point. For the location and height of the maximum in the susceptibilities 
one finds 

''max — iiOq , 
Xm,max — TTl^ — Ul^ ^ , (6) 

Xt,max = ct m-^* = ct -1)/^^ . 

For finite values of the quark mass the pseudo-critical couplings defined through peaks 
in Xm or xt can, of course, differ. The exponents Zt and Zm, which control the scaling 
behaviour of the thermal and chiral susceptibilities are also given in Table I. 



Similarly we obtain for the scaling behaviour of the chiral order parameter on the zero- 
mass critical point {t = 0) as well as on the line of pseudo-critical couplings {t = tmax), 

(«=c^my' • (7) 

In the case of spin models the analysis of cumulant ratios, which themselves are scaling 
functions, has proven to be useful. The simplest ratio one can consider in the case of QCD 
is the ratio of the second and first derivative of the free energy with respect to the quark 
mass, ie. the ratio 

As can be seen again from eq. (5) this ratio is a scaling function, which only depends on 
the combination y = th~^/^^ . A consequence of this is that A{t = 0, /i) is unique for 
all values of h (of course, modulo corrections from the regular part of the free energy). 
Outside the critical region the order parameter (ipip) is expected to depend linearly on the 
quark mass, {if^ij) = cq + ciruq with cq being non-zero for t < and zero otherwise. From 
this one obtains 

r 1 t <0 , h^O 
A{t, h) = Q{th-^/^^) = \l/5t = Q . (9) 

I t > , /i ^ 

Similarly one can analyze the scaling behaviour of the specific heat peak. Like in 
the case of the susceptibilities, the quark mass dependence of its location is given by tmax 
defined in eq. (6). However, unlike for the susceptibilities discussed above, the analytic 
calculations for the exponent a suggest that the specific heat does not diverge in the zero 
quark mass limit {a < 0) but scales as 

C'v^.max = Co + ci m~^" ; Za = a/pS . (10) 

The mass dependence thus is a subleading term. 



3) Determination of Critical Exponents on the Lattice 

In the following wc will discuss the observables from which information on the critical 
exponents a, f3 and d is extracted. These are various response functions which contribute 
to the specific heat and certain susceptibilities defined for two-fiavour QCD. In particular, 
we will study the behaviour of three susceptibilities - the chiral (Xm) and thermal (xt) 
susceptibilities, defined in eq. (1), as well as the Polyakov loop response function xl, 

XL = N'AiL') - (L)'} , (11) 



where L = J2x'^Yl.i^=i ^{xi) 6 denotes the average Polyakov loop on a lattice of 

size X Nt-. While the critical behaviour of the first two can be extracted from the 
structure of the free energy density, the latter is not directly related to it nor does it 
serve as an order parameter for a symmetry of the Lagrangian at finite quark mass. In 
fact, it is expected that (L) does not show any critical behaviour (divergence) in the limit 
of vanishing quark mass. Nonetheless, the existing Monte Carlo investigations show the 
presence of a pronounced peak in the Polyakov loop response function, which reflects the 
sudden onset of deconfinement in the chirally symmetric phase. 

On the lattice the chiral susceptibility Xm for flavours is obtained from matrix 
elements of the staggered fermion matrix, D = m^l + D^, as 



Xm — XO ~l~ Xconn 



(12) 



with 
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(13) 



Xconn — 

The first term, xoj gives the fiuctuations of the order parameter 

Nf 



(TrD" 



(14) 



The second term, Xconn, results from the connected part appearing in the second derivative 
of the fermion determinant. It is the integral over the connected part of a scalar propagator. 

The thermal susceptibility, xt, requires to calculate the derivative of the chiral con- 
densate with respect to the temperature, which in terms of the temporal lattice spacing, 
ttr, and the number of lattice points in that direction, A^,-, is given by T = (A^T-a,-)"^. 
Following the standard procedure for taking derivatives with respect to the temperature 
on the lattice [12] we obtain 



Xt 



1 d 



(15) 



where e denotes the energy density operator [12] 
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Here -Pct(t) = ^^a^^r^ Y2p Yli^p denotes the usual plaquette operator for space- 

like (timelike) plaquettes; 'ipDo'ip = jNfN~^ N~^TiDoD~^ while ^ = a^/ar is the ratio 
of spatial and temporal lattice spacing and the derivatives in eq. (16) should be evaluated 
for ^ = 1. For further details we refer to Ref. [12]. 

The calculation of xt gets complicated due to the fact that the derivatives of the 
couplings with respect to the temporal lattice spacing are only known perturbatively. These 
perturbative relations, however, are poor approximations to the exact derivatives in the 
coupling regime we are presently investigating. Moreover, for a complete determination of 
Xt one should also take into account the zero temperature part of the energy density, which 
still has to be subtracted in eq. (16). Fortunately, all this is not essential for the discussion 
of the scaling behaviour of xti ie. its divergence in the limit of vanishing quark mass. The 
zero temperature terms will not show any singular behaviour at the finite temperature 
phase transition point and the derivatives appearing in eq. (16) can be treated as constant 
factors. They will not influence the singular behaviour. We thus can study separately the 
scaling behaviour of the different terms contributing to Xt-, 

Xt,a^{^i^-Pa)-{^^){Pa) 

Xt,r = {^^-Pr)-{m{Pr) (17) 

Xtj = (V'V' • V'-DoV') - {ipij){ipDoip) 

Like in the case of the chiral susceptibility, the last response function, Xtj-, requires the 
calculation of contributions from connected diagrams. These have been omitted in the 
following analysis. Our results for Xm seem to justify this. Recall that, in contrast to Xm, 
we are here only interested in the divergent part of xt- For the same reason we also omit 
the contribution coming from the correlation of 'ipip with the mass term in eq. (16). Due 
to the additional mass factor this part will not contribute to the singular behaviour of Xt- 
Each of the other partial contributions to the thermal susceptibilities may diverge in the 
zero quark mass limit. In general, they even may diverge with different critical exponents, 

XM~m-^*- . (18) 

The critical exponent, Zt = {1 — j3)/ 135, which controls the divergence of Xt in the zero 
quark mass limit, is then given by the maximum of these three exponents 

Zt = maii{Zt,a: Zt,r: Ztj} . (19) 

Similarly we can investigate the scaling behaviour of the specific heat, defined in 
eq. (1). We can ignore all non-singular contributions and constant factors and simply 



investigate the scaling behaviour of the various response functions contributing to Cy; 



Here v denote a or r. We also left out the correlations of and ipDoip with the chiral 
condensate 'i/^i/j. These are given already in eq. (17) and, due to additional mass factors 
(see eq. (16)), certainly will not lead to a singular behaviour in the specific heat. From 
eq. (5) one expects that the quark mass dependence of the specific heat peak is controlled 
by the exponent Za = a/ (35. In the case of the 0{N) models in three dimensions the 
situation becomes, however, a bit more complicated due to the fact that the exponent a 
is expected to be negative [4] . In this case the specific heat is not expected to diverge, but 
rather should develop a cusp at Tg. Accordingly none of the response functions defined in 
eq. (20) is actually expected to diverge. They rather should approach a constant in the 
limit of vanishing quark mass, with the quark mass dependence described by eq. (10). 

4) Numerical Results 

We have performed calculations on an 8^ x 4 lattice, which may be considered to 
be quite a small lattice, yet, Ref. [5] shows that varying the spatial lattice extent from 
Ncr = 6 to 12 at A^T- = 4 surprisingly does not change the pseudo-critical temperature 
[8]. Likewise, increasing from 4 to 8 while keeping Nct/Nt between 2 and 4 also does 
not lead to apparent differences in the scaling behaviour of the pseudo-critical points. 
Nonetheless, it is important to get control over finite size effects. The present analysis 
should thus be considered as a first step towards a quantitative determination of critical 
exponents in QCD with light fermions. 

For the quark masses, we have chosen the values ruq = 0.075, 0.0375 and 0.02 in units 
of the lattice spacing a. These quark masses have been selected such that they can provide 
additional information on the quark mass dependence of pseudo-critical couplings beyond 
the existing data [2]. At each value of the quark mass we have performed simulations at 
several values of the gauge coupling, in the pseudo-critical region. The simulations 

have been carried out by means of the so-called hybrid R algorithm [13] . Guided by analyses 
of the hybrid Monte Carlo algorithm, we have scaled the step size according to the quark 
mass, dr = 0.05,0.04 and 0.03 (in the normalization of [13]) for = 0.075,0.0375 and 
0.02 respectively. The trajectory length was fixed to r = 1. For the conjugate gradient 
inversion we have settled at a required precision of ($ — DX)^/$^ < 0.5 x 10"-^^. The 




{ijDoij-P^)-{i^Doij){P^) 



(20) 



autocorrelation times have been determined to Tgxp < 100 trajectories right on the pseudo- 
critical couplings, the longest correlations being in the Polyakov loop, while away from the 
cross-over Tgxp is considerably smaller. At each value of the coupling we collected between 
5000 and 10000 trajectories which then have been analyzed using a reweighting with the 
density of states at intermediate values of the gauge coupling [14] . Errors were computed 
by the jack-knife procedure. 

4.1) Susceptibilities and the Specific Heat 

In Fig. 1 we show the Polyakov loop response function, xl-, for three values of the quark 
mass. Pronounced peaks are clearly visible. Although the location of the peaks clearly 
depends on m^, their heights do not have any significant quark mass dependence. This 
is consistent with the analysis presented in Ref. [8], which also suggests that the volume 
dependence of the peak height is small and distinctively different from the situation in the 
pure gauge theory. 

The behaviour of xl is also drastically different from that of the chiral susceptibility, 
Xmi which is shown in Fig. 2. The rise in the peak height with decreasing quark mass as well 
as the shift in the critical coupling is clearly seen. In Fig. 2 the two contributions, xo and 
Xconn are shown separately. The contribution from the connected part to the susceptibility 
is a slowly varying function in the critical region. The ratio Xconn/xo decreases with 
decreasing quark mass. However, even for ruq = 0.02 we find that Xconn contributes about 
30% to the value of Xm,max cind in the region of the zero quark mass critical temperature 
it gives the dominant contribution. This will be particularly important for our discussion 
of the chiral cumulant in section 4.3. 



Zm 


Zt,a 


Zt,T 




Zl 


0.79(4) 


0.63(7) 


0.63(7) 


0.65(7) 


0.05(6) 



Table II: Critical exponents controlling the divergence of the chiral (m), 
Polyakov loop (L) and various parts of the thermal susceptibilities in the limit 
of vanishing quark mass. Results are obtained from straight line fits to the data 
shown in Fig. 4. 

In Fig. 3 we show the partial contributions to the thermal susceptibility, defined in 
eq.(17). It is apparent that all three have a similar quark mass dependence. Note also 
that the peak position in all three cases coincides with the pseudo-critical couplings in xl 
or Xm- 



The peak heights of all susceptibilities are shown in Fig. 4 as a function of quark 
mass. A fit with a powerlike singular behaviour as indicated in eq. (18) yields the critical 
exponents given in Table II. While the exponent is compatible with being zero, we 
find that Zm is in remarkably good agreement with the 0{N) prediction. The agreement 
is not that good for the exponent of the thermal susceptibility, which comes out to be 
about 50% larger than expected for an 0{N) symmetric model. This, however, is not that 
unexpected as the exponent Zm is only sensitive to the magnetic exponent yh, while Zt is 
also related to yt- The latter gives directly the correlation length exponent yt = l/Su and 
thus is expected to be more sensitive to the spatial size of the system. We can rephrase 
our numerical results for Zm and Zt in terms of the exponents y^ and yt. This yields 

Vh = 0.83 ±0.03 , 

(21) 

yt = 0.69 ± 0.07 , 

which should be compared with the 0(4) values y^ = 0.828 and yt = 0.452. 

As has been discussed above the various response functions contributing to the specific 
heat are not expected to diverge in the zero quark mass limit. Results for three of the six 
different response functions defined in eq.(20) are shown in Fig. 5 for the three different 
quark mass values studied here. A comparison with the susceptibilities displayed in Figs. 3 
and 4 clearly shows that the increase in the peak height is considerably slower. However, we 
do not have any direct evidence for a finite limit of the peak height in the limit of vanishing 
quark mass. Results for the calculated peak values, Cy^^, and the corresponding pseudo- 
critical couplings are given in Table III. We have analyzed the scaling behaviour of these 
peaks in two ways, which clearly show the present uncertainties in the determination of 

(i) Assuming that the specific heat diverges in the zero quark mass limit, 

one can analyze the quark mass dependence of the specific heat response functions in 
the same way as the various components of the susceptibilities, ie. we fit the peak 
heights to eq. (10) with cq = 0. By construction this leads to Za > 0. We find, 
depending on the particular response function, values for Za between 0.28 and 0.34. 
The best fit result yields Za = 0.28(6). 

(ii) Assuming that the specific heat stays finite in the zero quark mass limit, 

one can fit the peak heights to eq. (21) with cq > 0. By construction this leads to 
Za < 0. As can be seen in Fig. 5 the strongest quark mass dependence is found in the 
purely gluonic response functions. They lead to Za = —0.05(2). 

We note that even by assuming a finite specific heat in the limit of vanishing quark 
mass our fits suggest a large value of the response functions in that limit. For instance we 



find C™^^(mq = 0) ~ 6. This indicates that these response function will still show quite 
a strong quark mass dependence also for quite small values of m^, which is related to the 
nearly vanishing magnitude of a (a = corresponds to a logarithmic singularity). The 
upper limit for a obtained in this analysis, a < 0.34, corresponds to an upper limit of 0.6 
for the thermal exponent yt- As can be seen from eq. (21) this is slightly lower than the 
result obtained from the analysis of the thermal susceptibility. 

Further calculations at smaller quark masses are thus particularly interesting for these 
response functions and are needed in order to extract the thermal exponent yt and to 
further clarify the nature of the scaling behaviour of the specific heat in the zero quark 
mass limit. 
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m 


Da 


Dt 


err 


aa 


TT 


0.075 


0.0535(53) 


0.0691(60) 


0.761(65) 


0.807(58) 


0.957(71) 




5.350(4) 


5.350(4) 


5.348(4) 


5.348(4) 


5.348(4) 


0.0375 


0.0722(83) 


0.0921(97) 


0.974(45) 


0.999(40) 


1.193(51) 




5.306(4) 


5.306(4) 


5.306(2) 


5.306(2) 


5.306(2) 


0.02 


0.0844(61) 


0.1040(66) 


1.165(62) 


1.170(53) 


1.409(70) 




5.282(2) 


5.282(2) 


5.282(2) 


5.282(2) 


5.282(2) 



Table III: The peak values of various response functions contributing to the 
specific heat and the corresponding values of the pseudo-critical couplings. 

4.2) Pseudo- Critical Couplings 

The locations of the peaks in the various susceptibilities or the specific heat, which 
have been discussed in the previous section, can be used to define pseudo-critical cou- 
plings. Although these pseudo-critical couplings may differ for different observables at 
non-vanishing values of the quark mass, we found that they agreed within our numeri- 
cal accuracy. We thus may quote a common pseudo-critical coupling for all the response 
functions analyzed by us. These are collected in Table IV. 

It is quite reassuring that the pseudo-critical couplings extracted from the location of 
the peak in Xm, the various components of xt and Cy agree this well. In particular this 
shows that the pseudo-critical couplings, extracted so far from the variation of the slope 
of {'ipip) as a function of (this corresponds to the location of the peak in xt,a(T)) can 



indeed be used to discuss the scaling behaviour of these couphngs as a function of the 
quark mass. 



Ill,, 




0.075 

0.0375 

0.02 


5.350(4) 
5.306(2) 
5.282(2) 



Table IV: Pseudo-critical couplings determined from the location of the peak 
in various response functions. Errors are obtained from a jack-knife analysis of 
the interpolation curves resulting from a reweighting with the density of states 
method. 

In Fig. 6 we show the quark mass dependence of the pseudo-critical couplings. Here 
we also include data from earlier simulations [6,7], which have been summarized in Ref. [2]. 
The pseudo-critical couplings have been fitted with the ansatz 

^ =co + cim^= , (22) 



where cq gives the critical coupling in the zero quark mass limit and Zc = 1/ PS. We find 

^r-- = 5.243 ±0.010 

''<-°^ . (23) 

— = 0.77 ±0.14 

Compared to the corresponding 0(4) value 1//35 = 0.55(2) and the 0(2) value 1/135 = 
0.60(1) the three parameter fit yields a larger value for this combination of exponents. 
However, also a two parameter fit with 1/(35 fixed to its 0(4) value still yields a good x^. 
We show this fit also in Fig. 6. The resulting zero quark mass critical coupling is found to 
be 6/(7^(0) = 5.222 ± 0.003. We also note that the result obtained for 1/135 is consistent 
with the determination of the thermal and magnetic exponents in the previous section and 
also suggests that in our present analysis we overestimate the value for yt. 

4.3) The Chiral Cumulant 

The most direct procedure to determine the critical exponent 5, which is independent 
of any ansatz for a fitting function, is given through an analysis of the chiral cumulant, 
defined in eq. (8). In the vicinity of the critical point this is a simple scaling function and 



curves for different quark masses cross in a unique point - the zero mass critical coupling, 
if subleading corrections from the non-singular part of thermodynamic quantities can be 
ignored. In Fig. 7 we show A(6/^^, m^) for our three different quark masses. For the largest 
mass the pseudo-critical region was quite far away from the zero quark mass critical region 
and we did not attempt to extend the calculations for that mass into this regime. The 
simulations for the two smaller quark masses, however, have been extended into this region 
and we see from Fig. 7 that they indeed cross at a coupling, which is consistent with the 
result obtained from our fit to the pseudo-critical couplings (Fig. 6). Moreover, we find 
that the value of A at the crossing point gives an astonishingly accurate estimate of 1/5. 
In fact, due to the rather weak dependence of A on the coupling Q/g^ in this region, we 
find that the exponent 1/5 is much better determined than the zero quark mass critical 
coupling. Taking for the latter the interval of values obtained from the two fits shown in 
Fig. 6. ie. 5.22 < Q/g1{Q) < 5.25, we get for the exponent 5, 

0.21 < 1/5 = 0.23 < 0.26 . (24) 

5) Conclusions 

We have studied various response functions which contribute to susceptibilities and 
the specific heat in two-fiavour QCD, simulated with staggered fermions on a lattice of size 
8^ X 4. We find that these observables are very well suited to localize the pseudo-critical 
couplings at non- vanishing values of the quark mass. At least at the moderate values of 
the quark mass analyzed here the height of the peaks in the susceptibilities is very well 
determined and can be used to investigate the scaling behaviour of these quantities. The 
analysis of the peak heights in the response functions as well as the quark mass dependence 
of the peak location yields independent observables, which can be used to determine the 
critical exponents a, j3 and 5. Additional information on 5 and the location of the zero 
quark mass critical point is obtained from the structure of the chiral cumulant. These 
exponents can be related to the two basic exponents yt and y/^, which characterize the 
scaling behaviour of the singular part of the free energy density. 

Our present analysis yields a magnetic exponent ym-, which is consistent with the 
value expected for a second order phase transition controlled by 0(4) exponents (eq. (21)). 
Our result for the thermal exponent, however, turns out to be about 50% larger than the 
corresponding 0(4) value. 

Clearly the errors are still quite large. However, the present results are consistent with 
our expectations based on the existence of a second order phase transition in the limit of 



vanishing quark mass. We find it particularly reassuring, that we obtain the weakest 
quark mass dependence in response functions which contribute to the specific heat. We 
also confirm a weaker quark mass dependence of the peak of the thermal susceptibility 
relative to that of the chiral susceptibility. Moreover, we find that the Polyakov loop 
response function is insensitive to changes in the quark mass. This suggests that, indeed, 
the chiral symmetry restoration is the driving mechanism for the finite temperature phase 
transition in two flavour QCD. 

We expect that in particular the thermal exponent yt is sensitive to finite lattice effects. 
The behaviour of the susceptibilities thus has to be studied for smaller quark masses and 
on larger lattices. 
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Figure Captions 

Figure 1: The Polyakov loop response function versus 6/(7^ for three values of the quark 
mass. Shown are results from a reweighting analysis with bins of length A (6/(7^) = 
0.002 close to the pseudo-critical coupling and 0.005 elsewhere. The couplings at 
which simulations have actually been performed are marked by filled symbols. 

Figure 2: The chiral susceptibility versus 6/(7^ for three values of the quark mass. Results 
from a reweighting analysis for the disconnected part xo cind the connected part 
Xconn defined in eq. (13) are shown separately. 

Figure 3: The three response functions which contribute to the thermal susceptibility de- 
fined in eqs. (15)-(17). 

Figure 4: Peak values of the Polyakov loop response function (xl), the chiral {Xm) and the 
thermal (xt) susceptibility versus the quark mass. 

Figure 5: Three response functions which contribute to the specific heat as defined in 
eq. (20). 

Figure 6: Pseudo-critical couplings, 6/(7^ (mg), versus h = NriTiq for A^,- = 4. The three 
new values analyzed here are shown as full diamonds. The other data points are 
taken from the collection given in Ref. [2] . The dashed curve is a three parameter 
fit with eq. (22). The solid curve is a two-parameter with this function and Zc 
fixed to the 0(4) value Zc = 0.55. 

Figure 7: The chiral cumulant, A, versus G/g"^ for three values of the quark mass. 
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